SI Methods 
Polymer model 

Model overview: We modeled chromatin as a fiber of monomers 
connected by harmonic bonds. Unless noted, each spherical monomer 
had a diameter of 15 nm and represented 500bp, or approximately 
three nucleosomes. A permanent loop was formed by connecting two 
monomers with a harmonic bond of the same strength as the bonds 
between all adjacent monomers. This permanently brings the loop 
bases into contact. Two such loops were formed in the two-loop 
simulations. A three-point interaction force was used to impose a 
bending energy and account for the rigidity of the fiber. To model 
volume interactions, monomers interacted via a shifted Lennard-Jones 
potential, which is a computationally fast purely-repulsive potential. 
Unless noted otherwise, the Lennard-Jones potential was truncated at 
U = 3kT as specified below to allow occasional fiber crossing. Simulated 
polymers were confined to a sphere at a given density and initialized 
from an unentangled polymer conformation. 

Polymer models were simulated with OpenMM - a high-performance 
GPU-assisted molecular dynamics software 

( https://simtk.org/home/openmm ). We used an in-house openmm- 
polymer library to efficiently set up polymer simulations with OpenMM, 
and to analyze simulation results, openmm-polymer is publicly 
available on the Bitbucket online repository: 
http://bitbucket.org/mirnylab/openmm-polymer 

Simulations were characterized by 4 parameters: loop size, number of 
loops, fiber stiffness, and system density. Total polymer length was 
always chosen to be approximately 33 * loop size. The initial 
conformation for all simulations was an unentangled polymer ring. 
Simulations for a phantom chain were performed by switching off 
inter-monomer Lennard-Jones interactions. Choice of parameters for 
various models is summarized in Supplemental Table Tl. 

Forces and Langevin Dynamics simulations 

Adjacent monomers were connected by harmonic bonds with a 
potential U = 25*(r - l) 2 (the energy is given in units of kT). The 
stiffness of the fiber was modeled by a three point interaction term, 
with the potential U = k*(l-cos(a)), where a is an angle between 
neighboring bonds, and k is a parameter controlling stiffness. A value 
of k=3 was used for most simulations; k=2 and k=4 were used for 
simulations with lower and higher stiffness; k=2 was used for 
simulations with the smallest (10 monomer) loop. 



Neighboring monomers interacted via a shifted Lennard-Jones (LJ) 
repulsive potential U = 4 * (1/r 12 - 1/r 6 ) + 1, for r<2 1/6 ; U=0 for r > 
2 1/6 (for details see (25)). To account for activity of type-II 
topoisomerase, we allowed fiber crossing by truncating the shifted LJ 
potential at an energy of E cu toff = 3 kT. For energy U more than 0.5 
Ecutoff, the LJ potential was modified as: Usoftened = 0.5 * E cu toff * (1 + 
tanh(2*U/E cu toff - 1))- To avoid numerical instabilities, the interaction 
radius r was truncated at r = 0.3 via: r tr uncated = (r 10 + 0.3 10 ) 1/10 , which 
introduced a negligible shift in the final softened potential. Spherical 
confinement was realized as a potential linearly increasing at a rate of 
k=5kT/mon when the radius was larger than the confinement radius. 

We simulated our model using Langevin Dynamics simulations, 
performing 80,000 blocks of 3000 time steps (240,000,000 time steps 
total). For the fiber lengths considered here, polymer simulations 
reached equilibrium in less than 1000 blocks; this was confirmed by 
observing that monomer displacement saturates at about 500 blocks. 
Polymer conformations starting with block 1000 were used for our 
simulations. We note that this study focuses on equilibrium aspects of 
chromatin loops and that simulated time is not specifically matched to 
the time-scale of enhancer-promoter interactions in vivo; a study of 
chromatin loop dynamics over time is beyond the scope of this study. 

An Andersen thermostat was used to keep the kinetic energy of the 
system from diverging. The time step was then chosen to ensure 
conservation of kinetic energy and lack of fiber crossing with the non- 
truncated Lennard-Jones potential. The absence of fiber crossing in 
this case was confirmed by the conservation of Alexander's polynomial 
for a 50000-monomer ring simulated at a high density of 0.85 for 
1,000,000,000 time steps. 

Initialization and starting conformations 

Since our simulations were performed at thermodynamic equilibrium, 
the starting conformation does not affect properties of the resulting 
heatmap; for simulations with fixed topology (i.e. no fiber crossing), 
only the topological state of the starting conformation is relevant. For 
simulations with or without fiber crossings, we initialized our 
simulations from an un-entangled polymer state created as described 
below. We started with a 4-monomer ring on a cubic lattice. We then 
chose one bond at random, and tried to extend the polymer at this 
location by two monomers, by making a bond into a kink. To do this, 
we considered another bond, obtained by shifting this bond by one in a 
random direction perpendicular to the bond (choosing one out of 4 
possible directions). If both locations of the shifted bond were free, the 



polymer was extended to incorporate this bond. For example, if a 
chosen bond was going in +z direction: ... -> (0,0,0) -> (0,0,1) -> ... , 
and we attempted to grow it in the -y direction (chosen randomly out 
of +x, -x, +y, -y), we would check positions (0,-1,0) and (0,-1,1). If 
both of them were free, the polymer sequence would be changed to ... 
-> (0,0,0) -> (0,-1,0) -> (0,-1,1) -> (0,0,1) -> ... . If at least one of 
them was occupied, selection of the random bond was repeated. The 
process was repeated till the polymer grew to the desired length. Since 
no polymer fibers can pass between the old bond and a kink, this 
process preserves the original topology and creates an un-entangled 
polymer. 

Calculating and analyzing heatmaps 

To obtain heatmaps, we first found all contacts within each polymer 
conformation. A contact was defined as two monomers being at a 
distance less than 2 monomer diameters. Contacts for all pairs of 
monomers were then put on a heatmap (i.e. a 2000-monomer 
polymer produced a 2000x2000 heatmap). When calculating contact 
frequency ratios for insulation and facilitation, averaging was 
performed over small regions of the heatmap to reduce sampling 
noise. Unless noted, we report the average value for insulation over a 
region of the heatmap, by averaging over monomers in the promoter- 
proximal third of the loop and over a +/-3 monomer enhancer- 
promoter separation. The range of insulation values for different 
positions in the promoter-proximal third of the loop is shown in Figure 
3B. For facilitation, we average over a region of the heatmap defined 
by a +/- 3 monomer enhancer-promoter distance and a +/- 
((enhancer-promoter genomic distance - loop length) / 6) monomer 
offset from a symmetric placement of an enhancer-promoter pair 
around the loop bases (e.g. +/- 6 monomers in Figure 2B). To 
calculate the contact frequency ratio, we used simulations without a 
loop to calculate the expected frequency; for all parameter values, two 
simulations without a loop and ten simulations with a loop were 
performed, with a newly generated starting conformation for each 
simulation. 





Supplemental Figure SI: Effects of chromatin fiber flexibility. 

(A) Heatmap on left displays log (total # of contacts) for simulations 
with a more flexible polymer and standard parameters: 30kb 
chromatin loop, 2% density, fiber crossing (topoisomerase activity). 
On the right is a heatmap for the less flexible polymer. In both cases, 
the loop features observed in Figure 2B are still present. (B) Bar plot 
shows insulation and facilitation: at the stiffness presented in the main 
figures, for a more flexible polymer, and for a less flexible polymer. 
(C) (top) shows a 20 monomer or lOkb stretch from a conformation of 
a more flexible polymer, (bottom) shows a 500 monomer or 250kb 
region from a conformation of a more flexible polymer. (D) Same as 
(C), but for a less flexible polymer. Note the smoother appearance of 
the chromatin fiber conformation. 
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Supplemental Figure S2: Effects of loop size. Schematics show 
insulation and facilitation arrangements including the enhancer 
(yellow), the promoter (blue), and the loop bases (orange) for the 
heatmaps below. In all cases the main qualitative features remain the 
same. (A) A 300kb by 300kb heatmap for a 30kb loop, as shown in 
Figure 2B. (B) A 300kb by 300kb heatmap for a smaller loop of length 
15kb. (C) A 300kb by 300kb heatmap for a larger loop of length 60kb. 
(D) A 50kb by 50kb heatmap for a very small loop of length 2.5kb with 
a more flexible fiber (k=2, see SI Methods). In this simulation only, 
each monomer represents 250bp rather than 500bp to account for 
local nucleosome destabilization and turnover. The 10 monomer loop 
considered here thus represents 2.5kb of a loosely-arranged chromatin 
fiber. Note the color of the map differs due to the smaller dynamic 
range in total number of interactions for this shorter chromatin fiber. 
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Supplemental Figure S3: Effects of phantom polymer chain. (A) 

Heatmap for phantom polymer chain with a 30kb loop, where 
insulation and facilitation arrangements are shown as in Figure 2B. The 
vertical and horizontal stripes of depleted interactions are almost non- 
existent, indicating dramatically reduced insulation. (B) Bar plot 
displays insulation and facilitation for the regular scenario (Figure 2B) 
on the left and the phantom chain on the right. Facilitation is slightly 
diminished, whereas insulation completely disappears. 
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Supplemental Figure S4: Effects of two consecutive loops. (A) 

Four sample polymer conformations from simulations of the two-loop 
system with loops (black) and loop bases (orange) highlighted. (B) 
Heatmap which shows log (total # of contacts) for the two-loop 
system. Each loop is 30kb in the 300kb by 300kb region shown. The 
four red dots closer to the diagonal are the direct interaction of the 
loop bases from the formation of two loops. The two, weaker, red dots 
further from the diagonal are the interaction between the base at the 
start of the first loop and at the end of the second loop. The horizontal 
and vertical stripes of darker blue are indicative of strong insulation. 
Annotations show two loops formed from three bases (orange) along 
with the insulation and facilitation enhancer-promoter placements. (C) 
Schematics of enhancer-promoter arrangement for the two-loop 
system (top) insulation, indicated the red "X", (bottom) facilitation. 



Name 


Figure # 


Fiber crossing 


Density 


Flexibility (k) 


Loop size (kb) 


Excluded volume 


# of loops 


1 monomer = x bps 


Standard 


2,3A,3B 


yes 


0.02 


3 


30 


yes 


1 


500 


Without topo-ll 


3C 


no 


0.02 


3 


30 


yes 


1 


500 


Low density 


3D 


yes 


0.01 


3 


30 


yes 


1 


500 


Density 0.05 


3D 


yes 


0.05 


3 


30 


yes 


1 


500 








0.1 


3 


30 


yes 


1 




High density 


3D 


yes 


0.2 


3 


30 


yes 


1 


500 


More flexible 


SI 


yes 


0.02 


2 


30 


yes 


1 




Less flexible 


SI 


yes 


0.02 


4 


30 


yes 


1 


500 


Smaller loop 


S2 


yes 


0.02 


3 


15 


yes 


1 




Larger loop 


S2 


yes 


0.02 


3 


60 


yes 


1 


500 


Smallest loop 


S2 


yes 


0.02 


2 


2.5 


yes 


1 




Phanton chain 


S3 


NA 


0.02 


3 


30 


no 


1 


500 



Supplemental Table Tl: List of parameter values for all presented 
simulations. 



Supplemental Movie Ml: Langevin Dynamics of 30kb permanent 
loop formed in 1Mb regions of chromatin fiber, as in Figure 2. The 
polymer is colored according to the facilitation arrangement, where the 
loop (black) occurs in the region between the enhancer (yellow) and 
the promoter (blue). The movie is presented at a rate of 1000 
simulated time-steps per one second of real time; every seventh frame 
of the movie corresponds to a computationally obtained conformation, 
with quadratic interpolation performed between subsequent 
conformations (done to avoid high-frequency fluctuations between 
neighboring frames). The total simulation time for each run of each 
parameter set was approximately 5,000 times longer than the part of 
simulation displayed in this movie. 



